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Abstract: Voltage stability is a dynamic phenomenon and in order to tackle the voltage stability problem more realistically, 
dynamic representation of all power system components is necessary. The influence of load on dynamic stability is 
significant. In most studies, load representation of power systems are considered as constant power type but in real life there 
are different types of loads as industrial, agricultural and residential etc. Loads in a power system have to be represented 
realistically by adding static and dynamic types of loads to the system. It is known that a major portion of load is induction 
motors. In this paper, microgrid test system will be expended to a point that various combinations of static and dynamic 
loads can be considered as loads. After modeling part, effectiveness of the microgrid voltage stabilizer (MGVS) developed 
for microgrids that have only constant power type of loads will be studied under various combinations of static and dynamic 
loads. 
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I. NOMENCLATURE 

Subscript m - imaginary axis; 
Subscript r — real axis; 

Vo> Po > Qo ~ initial conditions of the system, (Voltage, active power, reactive power respectively); 
% to a 6 — Parameters of ZIP load model; 
K p f, K q f — Frequency sensitivity parameters; 

e —First cage voltage; 

H b — Base radian electrical frequency; 

a- Slip; 

x R — Cage rotor reactance; 

r R — Cage rotor resistance; 

T M — Mechanical torque applied at the shaft; 

T e — Electrical torque; 

x m — Magnetization reactance; 

x s — Stator reactance; 

r s — Stator resistance; 

H — Shaft inertia constant; 

V t des ~ Desired voltage 

Vi dyn ~ Dynamic voltage 

AVi err — Per unit (pu) difference between the desired voltage and the dynamic voltage 

AV err — Total voltage deficiency 

K — Gain constant 

T lt T 2 — Time constants 

Vmgvs ~ Output of the controller 

oq, oc 2 , ... . ocj —The weighting factors for load buses 

II. INTRODUCTION 

Load representation is a critical issue on approaching power system modeling realistically. Loads are still 
considered as one of the most uncertain components of a power system to model because of their randomness, different 
timescale properties and their statistical nature. The results of the stability studies of power systems depend on how the load 
models represent the real load types. Several studies have shown the crucial effect of load representation in voltage stability 
studies [1] - [2]. Induction motors (IM) are used for representing the dynamic part of the load models. Depending on the type 
of the user profile, between 50% and 70% of the entire load consists of three phase IMs [3]. Static Polynomial (ZIP) load 
types are also commonly used and can be classified into constant impedance (Z), constant current (I) and constant power (P) 
load. The power has a quadratic dependence on voltage for Z load. For an I load, this dependency is linear, and power is 
independent of voltage changes for P type of load. 

In addition to load representation and modeling, increasing power demand is another leading issue exists for utility 
companies. More demand on power is stressing the generation system capabilities and transmission lines. Low efficiency of 
central plants and transmission and generation losses plus frequent power outages cost the United States hundreds of billion 
dollars per year [4]. Researchers are looking for alternatives that can fix those pricey problems without adding new 
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transmission lines. Integration of the distributed energy sources to main grid can be a possible answer for these questions but 
this concept has its own shortcomings. Therefore, researchers and scientists proposed a network called "Microgrid". 
Several blackouts have been associated with voltage stability problems in a power system [5] -[6]. A fast voltage collapse can 
be avoided by using available dynamic reactive power capabilities sufficiently. Transferring the reactive power within long 
distances causes massive voltage drop for typical power systems. But this is not the case for microgrid systems. Sustaining 
the dynamic voltage stability for microgrids can be possible with coordinated compensation of reactive power sources 
because electrical distance between the loads and the sources of the reactive power are relatively short. This concept has 
been used for designing Microgrid Voltage Stabilizer (MGVS) [7]. 

In this study, a 21- bus microgrid system, as shown in Fig. 1, run by diesel engine generators (DEG) will be used 
and load types will be considered as IM and ZIP loads to make the microgrid environment as realistic as possible to 
investigate voltage stability issues. After adding the IM and ZIP loads to microgrid model, effectiveness of the MGVS will 
be investigated for all load types. 

III. MICROGRID 

A microgrid is a distributed energy system which is a part of a large power system and it is supported by one or 
more distributed generation (DG) units. Microgrids are almost 85% efficient and have combine heat and energy applications 
and lesser transmission losses [8]. Its capability of operating in parallel with the grid or being able to operate in islanding 
mode during power outages and disturbances provides higher flexibility and reliability of operation. In islanding mode, 
microgrids retain power availability, avoid lost productivity and blackouts. 

The aim of a microgrid is to provide a value to both customer and utility by supplying to local loads [9]. Higher 
power quality, reduction in environmental pollutants, higher reliability of power distribution, and decreasing power line 
congestion can be listed as some of the benefits of the integration of microgrids into utility grid. 

IV. MODELING OF LOADS 

Previous studies on modeling of the microgrid itself [10] is extended to a point that all of its loads at the microgrid 
are converted to static and dynamic type of loads from constant power loads. In this section, modeling of dynamic and static 
loads will be presented. 



1. STATIC POLYNOMIAL (ZIP) LOAD 

ZIP Load Model is one of the oldest load representations. It is called ZIP because it is a combination of Z, I and P 
type of loads. Following equations represent the static polynomial load model [11] 



P=P 0 [a 1 g +a 2 g+a 3 

r 2 

Q= Q ° hS +as S +a6 



a + K pf An 

(l + K qf Af) 



(1) 
(2) 



2. INDUCTION MOTOR 

In this study, third order single cage induction motor has been used [12]. The single-cage induction motor's 
simplified electrical circuit can be seen in Fig. 2. The equations are formulated in terms of real and imaginary axes, with 
respect to the network reference angle. 
The network and stator machine voltage is: 
v r = —V sin 6 (3) 



v m = V cos 6 



The power absorptions are: 



(4) 
(5) 
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Figure 1. 21- Bus Microgrid System. 
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Figure 2. Simplified electrical circuit for single cage IM. 



In terms of the voltage behind the stator resistance r s , the differential equations are [12] 

d ' n ' e r +(xo~ x )i ? 

— ^ r lii ) (J6 m 

d 

dt ^ 



= —(2 b ae r — 



e m +(x 0 -x )i, 



(7) 
(8) 



The link between state variables, currents and voltages is: 

d 

(10) 



X^-y G y* ^5 ^7* I ^777, 

d 

^5 X l r 

where x 0 , x and T 0 can be obtain from the motor parameters. 

*0 = x s + (11) 
x Rl x m (12) 

(13) 



x = x s + 



7"o = 



XRl+Xm 



n b r R1 

At last, the mechanical equation is: 

d T m (o)-T e 
— (J = 



dt 



2// 



m 



where the electrical torque is: 
7g 6? r i r ~h 6? m i m 



(14) 



(15) 



V. MICROGRID VOLTAGE STABILIZER 

The MGVS has a similar functionality compared to a power system stabilizer (PSS) in terms of approaching voltage 
stability problem of a power system. The PSS gives an input to the excitation system of a generator to bring voltage stability 
to a power system. In order to prevent any voltage collapse, The MGVS gives an input to reactive power loops of DGs or the 
excitation systems, which lets DGs to kick in more reactive power into the microgrid. By implementing this method, using of 
costly dynamic reactive sources like capacitor banks, SVC or STATCOM can be avoided. The MGVS model and its 
simplified version can be seen in Fig. 3 and Fig. 4 respectively. 

The MGVS input is a measurement of the per unit (pu) difference (AVJ err ) between the desired voltage (VJ d es ) an d 
the dynamic voltage (Vj d yn )• This voltage deficiency is calculated for all the load buses [7]. 

AV . ^ = V ides-V idyn . = ^ _ _ l (lg) 

V i des 

Based on the importance of the bus, weighting factors for all buses are defined. In order to get a total voltage 
deficiency (AV err ) of the system, a weighted average of AV t err is taken. A lead/lag block consisting of gain constant (K) and 
time constant 7\ and T 2 takes the AV err as an input. As shown in Fig. 3, V MGVS is the output of this MGVS controller [7]. 

a T7 oCjAVi err +°c 2 AV 2 err + ■■■■«/AV r z 

err 

L\ V err 



(17) 



oq+ oc 2 H hoc/ 

The weighting factors for load buses are oq, oc 2 , ... . oq. For this study, all oc values are equal to each other because 
loads at all buses are considered equally important. 

The block diagram in Fig. 4 can be implemented in MATLAB. The corresponding differential equations 
representing the MGVS are given below. 

(18) 
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Figure 3. Microgrid Voltage Stabilizer 
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Figure 4. Simplified Microgrid Voltage Stabilizer Model 

VI. SIMULATION AND RESULTS 

After modeling of the microgrid system and its components, voltage stability of the system will be studied by 
implementing a three phase short circuit fault to a selected bus. In this paper, MATLAB programming environment is used 
to model the 21 -bus microgrid system that has three diesel engine generators and 6 loads and to implement the fault 
conditions [13]. Fig. 1 shows the details of the microgrid system used for this study. There will be two case studies. In first 
case study, load types will be considered as either 100% ZIP or 100% IM type and response of the system after disturbances 
will be evaluated. For the second case study, load types will be considered as 50% ZIP and 50% IM type. The results will be 
used to compare the cases of with and without presence of the MGVS. This will help us to understand the behavior of 
MGVS with ZIP and IM type of loads and its effectiveness. 



Case study 1: Three phase short circuit fault at Bus 7 with 100% ZIP and 100% IM Loads 

In this case study, a three phase short circuit fault is applied to Bus 7 for 1.5 sec. The disturbance starts at 0.5 sec. 
and overall simulation time is 7 sec. Bus Voltages, Active Power and Reactive Power for Bus 15 have been plotted 
separately for various load types with and without presence of MGVS as listed below. 

1) 100% I M 4) 100% P with MGVS control 

2) 100% I M with MGVS control 5) 100% Z 

3) 100% P 6) 100% Z with MGVS control 

In each case, only one type of load is considered and all results are shown in one plot for comparison. Aim for 
having only one type of load is for observing the system behavior for those loads individually. Depending on load types, the 
voltage at Bus 15 drops at the time of fault to a point that is all different for each case. As it can be seen from the Fig. 5, the 
biggest drop took place with 100% P load and least drop was observed with 100% Z load. Without presence of MGVS 
voltage drops around 0.92 pu and stay around that point. The MGVS works effectively for all load types to compensate the 
voltage drop. Voltage turns back to nominal point in 1.5 s. The results show that when the MGVS is active, the response of 
system after the fault for covering the voltage is alike for all load types. But when we look at Fig. 6 for active power and Fig. 
7 for reactive power response of the system, they are all different for each case. 

Z and P loads are the upper and lower limit for voltage dependency of active and reactive power of ZIP load 
respectively. Any other percentage of these Z-I-P load types will be in that range. Active and reactive power at all load buses 
for 100% P load is constant because they are independent from voltage changes. Largest decrease on producing active and 
reactive power happened for 100% Z load because its dependence of voltage is quadratic. At the times of fault and the fault 
clearing, 100% IM load has instant pick values for reactive power as in Fig 7. With presence of the MGVS, there is a 
significant increase for reactive power and active power for all load types. It is normal to observe this behavior because 
MGVS is improving the voltage profile of the system by using the reactive power compensation. 

Case study 2: Three phase short circuit fault at Bus 7 with 50% ZIP and 50% IM Loads 

Depending on the type of the user profile, between 50% and 70% of the entire load consists of three phase IMs. For 
this case study, 50% ZIP and 50% IM Loads is considered in order to make the load profile more realistic and suitable for 
real life applications. A Three phase short circuit fault is applied to Bus 7 for 1.5 s. Total simulation time is 7 s. Bus 
Voltages, Active Power and Reactive Power for Bus 15 have been plotted for various load types for the case of with and 
without presence of MGVS as listed below. Results are shown in one plot for comparison for each case. 



1) 50% P + 50% IM with MGVS control 

2) 50%P + 50%IM 



3) 50% Z + 50% IM with MGVS control 

4) 50%Z + 50%IM 
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Figure 5. Comparison of Bus Voltage at Bus 15 for Case Figure 6. Comparison of Active Power at Bus 15 for Case 

Study 1. Study 1. 
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Figure 7. Comparison of Reactive Power at Bus 15 for 

Case Study 1 . 



Figure 8. Comparison of Bus Voltage at Bus 15 for Case 

Study 2. 



In this case study, we have 50% fixed IM load. 50% ZIP load considered as 50% Z and 50% P type of load in order 
to see the range of possible variations. Depending on Z and P loads, the voltage at Bus 15 drops to a point that is slightly 
different for each case. As it can be seen from Fig. 8, voltage drop is more for the case having P load compared to the case 
having Z load. Without presence of MGVS voltage drops around 0.92 pu and stay around that point. When the MGVS is 
active, voltage turns back to nominal point in 1.5 s. and there is an 8% improvement in bus voltages. But when we look at 
Fig. 9 for active power and Fig. 10 for reactive power response of the system, they are all different. As it was observed at the 
case study 1 , largest drop for active and reactive power occurred at the case having Z load because its dependence to voltage 
is quadratic. At the times of fault and fault clearing, instant pick values for reactive power can be seen because of IM load 
for all cases. With presence of the MGVS, there is a significant increase for reactive power and active power for all load 
types. MGVS is using the available reactive power capabilities to improve the voltage profile of the system. 



VII. CONCLUSION 

Modeling of static and dynamic loads of a 21 -bus microgrid test system and successful implementation of three 
phase short circuit fault with variable ZIP load and dynamic IM load has been presented using MATLAB. By using this 
simulation study, it was observed that the static ZIP and IM load models can have different load power characteristics 
depending on the relation with the voltage. Most of the real load types that are being tried to model have all variety of load 
types like industrial, residential and agricultural. Combining ZIP and IM loads made microgrid system more realistic because 
this is the common practice in real life. By changing the percentage of the ZIP load parameters and overall percentage of IM 
load, desired load characteristics can be accomplished. By using this test system, effectiveness of MGVS on voltage stability 
enhancement has been investigated. After testing MGVS with pure ZIP, pure IM and combination of ZIP+IM loads, results 
showed that MGVS is kicking enough reactive power to the system to prevent voltage collapse and it is enhancing the 
voltage profile of load buses for each type of loads successfully. 
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Figure 9. Comparison of Active Power at Bus 15 for Case 

Study 2. 
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